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Abstract 

Background: The inflammatory bowel diseases (IBD), Crohn's disease (CD) and ulcerative colitis (UC), result from the 
combined effects of susceptibility genes and environmental factors. Polymorphisms in genes regulating inflammation may 
explain part of the genetic heritage. 

Methods: Using a candidate gene approach, 39 mainly functional single nucleotide polymorphisms (SNPs) in 26 genes 
regulating inflammation were assessed in a clinical homogeneous group of severely diseased patients consisting of 624 
patients with CD, 411 patients with UC and 795 controls. The results were analysed using logistic regression. 

Results: Sixteen polymorphisms in 13 genes involved in regulation of inflammation were associated with risk of CD and/or 
UC (p<0.05). The polymorphisms TLR2 (rs1816702), NFKB1 (rs28362491), TNFRSF1A (rs4149570), IL6R (rs4537545), IL23R 
(rsl 1209026) and PTPN22 (rs2476601) were associated with risk of CD and the polymorphisms TLR2 (rs1816702), TLR4 
(rs1554973 and rs12377632), TLR9 (rs352139), LY96 (rsl 1465996), NFKBIA (rs696), TNFA (rs1800629), TNFRSF1A (rs4149570), 
IL10 (rs3024505), IL23R (rsl 1209026), PTPN22 (rs2476601) and PPARG (rs1801282) were associated with risk of UC. When 
including all patients (IBD) the polymorphisms TLR2 (rs4696480 and rsl 81 6702), TLR4 (rsl 554973 and rsl 2377632), TLR9 
(rs187084), TNFRSF1A (rs4149570), IL6R (rs4537545), IL10 (rs3024505), IL23R (rsl 1209026) and PTPN22 (rs2476601) were 
associated with risk. After Bonferroni correction for multiple testing, both the homozygous and the heterozygous variant 
genotypes of IL23R G>A(rs1 1209026) (OR CD ,adj: 0.38, 95% CI: 0.21-0.67, p = 0.03; OR| BD , a dj 0.43, 95% CI: 0.28-0.67, p = 0.007) 
and PTPN22 1858 G>A(rs2476601) (OR CD ,u'nadj 0.54, 95% CI: 0.41-0.72, p = 7*10" 4 ; OR| BD , unadj : 0.61, 95% CI: 0.48-0.77, 
p = 0.001) were associated with reduced risk of CD. 

Conclusion: The biological effects of the studied polymorphisms suggest that genetically determined high inflammatory 
response was associated with increased risk of CD. The many SNPs found in TLRs suggest that the host microbial 
composition or environmental factors in the gut are involved in risk of IBD in genetically susceptible individuals. 
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Introduction 

Chronic inflammatory bowel diseases (IBDs), Crohn's disease 
(CD) and ulcerative colitis (UC), are complex diseases that result 
from the interaction of numerous genetic and environmental 
factors [1]. 

Genetic association studies have identified innate immunity as a 
critical component in the development of IBD. Until now, more 
than 163 IBD susceptibility polymorphisms have been confirmed, 
most of which are associated with both CD and UC, by candidate 
and genome wide association studies (GWAS) [2—11]. However, 
these polymorphisms have been estimated to only account for 20% 
of the genetic heritage involved in IBD [12]. 

Functional polymorphisms in genes in the inflammatory 
pathways may explain some of the genetic heritage involved in 
IBD. The transcription factor NFkB is a central regulator of 
inflammation. NFkB can be activated by Toll like receptors 
(TLRs). TLRs recognize pathogen-associated molecular patterns 
(PAMPs) that are broadly shared by pathogens but distinguishable 
from host molecules such as bacterial or viral DNA, flagellin or 
lipopolysaccharide (LPS). The TLRs initiates a kinase cascade that 
ultimately activates the IKK-complex, which phosphorylates and 
degrades the NFkB inhibitor IrcBa. NFkB is shuttled from the 
cytosol to the nucleus where it initiates expression of pro- and anti- 
inflammatory cytokines including TNF-a, IL-6 and IL-10 [13]. 

To identify susceptibility loci we assessed 39 mainly functional 
polymorphisms in genes involved in inflammation, particular in 
the NFkB pathway, in a homogeneous Danish cohort of 624 
patients with severe CD, 411 patients with severe UC and 795 
healthy controls. The candidate gene approach using functional 
polymorphisms allows interpretation of the underlying biological 
mechanisms based on increased or decreased gene expression or 
protein activity. 

Functional polymorphisms in genes in the inflammatory 
pathways involved in regulation of the NFkB pathway (TLR2, 
TLR4, TLR5, TLR9, LT96, CD14, MAP3K14, SUM04, NFKBIA 
and NFKB1), TNF-ot signaling (TJVFA, TJVFRSF1A and TMFAIP3), 
cytokines regulated by NFkB (IL1B, IL1RN, IL6, ILK), IL17A and 
IFJVG) and other genes involved in regulation of inflammation 
[IL4R, IL6R, IL23R, TGFB1, PTPN22, PPARG and NLRP3) were 
studied, as genetically determined variation in the inflammatory 
pathways may be associated with severe disease among patients 
with CD and UC. 

Materials and Methods 

Cohort 

A prior anti-TNF naive Danish cohort of patients with IBD was 
established. In short, blood samples retrieved as part of the routine 
screening for latent Mycobacterium tuberculosis at Statens Serum 
Institut (SSI, Copenhagen, Denmark) and the Department of 
Respiratory Diseases B or the Department Clinical Microbiology, 
Aarhus University Hospital (Aarhus, Denmark) were collected 
from 01.09.2009 to 30.03.2011 (9217 patients). Patients with 
intestinal diseases (ICD-10 code K50-K63) were identified by 



linking the unique personal identification number of Danish 
citizens (CPR-number) from each blood sample with the National 
Patient Registry (2659 patients). Patient records from 18 medical 
departments were examined (1378 patients) and identified 1035 
ethnic Danish patients with IBD where blood and clinical data 
were available. The patients either received or were considered 
candidates to anti-tumor necrosis factor-a (TNF-a) therapy 
(infliximab or adalimumab). The control group consisted of 795 
healthy blood donors recruited from Viborg, Denmark [4] . 

Selection of polymorphisms 

A candidate gene approach was used with focus on polymor- 
phisms in the TNF-a and NFkB pathways. In addition, 
polymorphisms in genes which have been shown to be associated 
with CD and/or UC, polymorphisms in inflammatory cytokines 
and polymorphisms in TLR2 and TLR4 were included [14]. 

Functional polymorphisms in relevant genes were found by 
searching pubmed with "polymorphism AND gene-name AND 
(reporter gene OR luciferase OR ELISA OR enzyme-linked 
immunosorbent assay OR RT-PCR OR reverse transcriptase 
PCR OR EMSA OR electrophoretic mobility shift assay OR flow 
cytometry)". 

Genotyping 

For patients with IBD the DNA was extracted from cryopre- 
served blood clots by using the Maxwell 16 Blood purification kit 
(Promega) according to the manufacturers' instructions with a 
median yield of 4.90 ug (range 0.8-25 ug) pr 300 ul total blood 
[15]. For the healthy controls, DNA was extracted from EDTA- 
stabilized peripheral blood by either PureGene (Qiagen, Hilden, 
Germany) or Wizard Genomic (Promega, Madison, Wisconsin, 
USA) DNA purification kit according to the manufacturers' 
instructions [4]. Competitive Allele-Specific Polymerase chain 
reaction (KASP), an end-point PCR technology, was used by LGC 
Genomics for genotyping (LGC Genomics , Hoddesdon, United 
Kingdom) (http://www.lgcgenomics.com/). The SNPs studied 
were TLR2 (rs4696480, rsl816702, rsl 1938228, rs3804099), 
TLR4 (rsl2377632, rs5030728, rsl554973), TLRS (rs5744168), 
TLR9 (rsl87084, rs352139), LT96 (MD-2) (rsl 1465996), CD14 
(rs2569190), MAP3K14 (NIK) (rs7222094), SUM04 (rs237025), 
NFKBIA (IkBoc) (rs696, rsl7103265), NFKB1 (NFkBI) 
(rs28362491), TJVFA (TNF-a) (rsl800629, rsl800630, rsl799724, 
rs361525), TMFRSF1A (TNFR1) (rs4149570), TJVFAIP3 (A20) 
(rs6927172), IL1B (IL- 1 P) (rsl 143623, rs4848306, rsl 143627), 
IL-1RN (IL-1RA) (rs4251961), IL4R (rsl805010), IL6 
(rsl0499563), IL6R (rs4537545), IL10 (rsl800872, rs3024505), 
IL17A (rs2275913), IL23R (rsl 1209026), IFNG (IFN-y) 
(rs2430561), TGFB1 (TGF-pJl) (rsl800469), PTPM22 (rs2476601), 
PPARG (PPAR-y) (rsl 80 1282) and NLRP3 (rs46 12666). 

Genotyping of TNFA (TNF-a) -857 C>T (rsl799724) and - 
863 OA (rsl800630) failed due to their close proximity to each 
other. All genotyping of -857 C>T (rsl799724) either failed or 
were erroneously genotyped as homozygous wild type when the 
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patients were carriers of the AA genotype of —863 OA 
(rs 1800630) due to genotyping bias. 

The 39 genotypes were replicated in 94 randomly selected 
samples and yielded >99% identical genotypes. 

Statistical analysis 

Logistic regression was used to compare genotype distributions 
among patients with CD, UC and IBD versus healthy controls 
(Table SI and S2). Crude odds ratio and odds ratio adjusted for 
age, gender and smoking status were assessed. A chi-square test 
was used to test for deviation from Hardy-Weinberg equilibrium in 
the healthy controls and for haplotype analysis (Table S3, S4, S5). 

Statistical analyses were performed using STATA version 1 1 
(STATA Corp., Texas, USA). 

Ethics statement 

The study was conducted in accordance with the Declaration of 
Helsinki and was approved by the Regional Ethics Committees of 
Central (M20100153) and Southern (S-20120113) Denmark and 
the Danish Data Protection Agency of Central (RM: J. 2010-41- 
4719) and Southern (RSD: 2008-58-035) Denmark. The Ethics 
Committees gave suspension for obtaining written informed 
consent. 

Results 

Study population 

Characteristics of the Danish patients with CD, UC and healthy 
controls are shown in Table 1. 



The genotype distributions among the healthy controls deviated 
from Hardy-Weinberg equilibrium for TLR2 (—16934 A>T 
(rs4696480)) (p = 0.02), TLR4 (rsl554973 T>C) (p = 0.03), TLR9 
(1174 OA (rs352139)) (p = 0.02), and TGFB1 (-509 OT 
(rs 1800469)) (p = 0.02). None of the deviations remained statisti- 
cally significant after correction for multiple testing. 

Polymorphisms associated with risk of CD 

The homozygous variant genotype of TLR2 OT (rsl816702) 
(OR adj : 2.80, 95% CI: 1.03-7.62, p = 0.04), TNFRSF1A -609 O 
T (rs4149570) (OR ildj : 1.84, 95% CI: 1.19-2.84, p = 0.01) and 
IL6R OT (rs4537545) (OR adj : 1.73, 95% CI: 1.12-2.66, 
p = 0.01) were associated with increased risk of CD. Both the 
homozygous and the heterozygous variant genotypes of NFKB1 — 
94ins/del (rs28362491) (OR unadj : 0.80, 95% CI: 0.65-1.00, 
p = 0.05), IL23R OA (rsl 1209026) (OR ildj : 0.38, 95% CI: 
0.21-0.67, p = 9*10" 4 ) and PTPM22 1858 OA (rs2476601) 
(OR adj : 0.57, 95% CI: 0.39-0.83, p = 4*10" 3 ) were associated with 
reduced risk of CD (Table SI and S2). 

After Bonferroni correction for multiple testing both the 
homozygous and the heterozygous variant genotypes of IL23R 
OA (rsl 1209026) (OR adj : 0.38, 95% CI: 0.21-0.67, p = 0.03) 
and PTP.N22 1858 OA (rs2476601) (OR unadj : 0.54, 95% CI: 
0.41-0.72, p = 7*10 4 ) were associated with reduced risk of CD. 

The variant allele of the polymorphisms have been shown to 
increase TLR2 levels {TLR2 OT (rsl816702)), increase 
TNFRSF1A expression (TNFRSF1A -609 OT (rs4149570)), 
increase IL-6r and IL-6 levels (IL6R OT (rs4537545)), decrease 
NF-kB p50 subunit expression (NFKB1 -94ins/del (rs28362491)), 



Table 1. Description of the study participants. 







Crohns Disease (CD) 


Ulcerative Colitis (UC) 


Controls 




(n = 624) 


(n = 411) 


(n = 795) 


Gender: n (%) 


Male 


272 (44) 


201 (49) 


411 (52) 


Female 


352 (56) 


210 (51) 


384 (48) 


Age: 


Median (5%-95%) 


37 (20-67) 


42 (20-72) 


43 (23-60) 


Age at diagnosis: 


Median (5%-95%) 


25 (14-59) 


33 (15-67) 




Smoking habits: n (%) 


Smokers 


1 78 (29) 


30 (7) 


207 (26) 


Former smokers 


64 (10) 


86 (21) 


392 (49) 


Never smokers 


156 (25) 


102 (25) 


189 (24) 


Data not available 


226 (36) 


193 (47) 


7(1) 


Location UC: n (%) 


Proctitis (El) 




53 (13) 




Left side (E2) 




1 83 (45) 




Extensive (E3) 




134 (33) 




Data not available 




41 (10) 




Location CD: n (%) 


Colonic (L2) 


208 (33) 






Ileal (LI) 


1 72 (28) 






lleocolonic (L3) 


210 (34) 






Data not available 


34 (5) 







doi:1 0.1 371 /journal.pone.009881 5.t001 
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decrease IL-17 serum levels {IL23R OA (rsl 1209026)) and 
decrease TNF-ot serum levels (PTPM22 1858 G>A (rs2476601)), 
respectively (Table 2). 

Thus, polymorphisms associated with higher TLR2 levels 
(TLR2 C>T (rsl816702)), increased TMFRSF1A expression 
(TMFRSF1A -609 G>T (rs4149570)) and higher IL-6r and IL-6 
levels (IL6R OT (rs4537545)) were associated with increased risk 
of CD. In addition, decreased IL-17 serum levels (IL23R G>A 
(rsl 1209026)) and decreased TNF-ot serum levels (PTPM22 1858 
G>A (rs2476601)) were associated with reduced risk of CD. The 
NF-kB p50 subunit can act both pro-inflammatory as part of the 
p50/p65 complex or anti- inflammatory as p50 homodimer 
(Table 2). A lowered NF-kB p50 subunit expression (JVFKB1 - 
94ins/del (rs28362491)) could be a differential lowering of the 
anti-inflammatory response which was associated with reduced 
risk of CD. 

Polymorphisms associated with risk of UC 

The homozygous variant genotype of TMFRSF1A —609 G>T 
(rs4149570) (OR adj : 1.85, 95% CI: 1.16-2.97, p = 0.01) and 
PPARG C>G (rsl801282) (OR unadj : 2.12, 95% CI: 1.01-4.45, 
p = 0.05) and both the homozygous and the heterozygous variant 
genotypes of TZJtf OT (rsl816702) (OR adj : 1.63,95% CI: 1.13- 
2.36, p = 0.01), TLR4 T>C (rsl2377632) (OR adj : 1.42, 95% CI: 
1.01-2.00, p = 0.04), MFKBIA 2758 A>G (rs696) (OR adj : 1.45, 
95% CI: 1.03-2.05, p = 0.03) and ILK) OT (rs3024505) 
(OR unadj : 1.42, 95% CI: 1.10-1.82, p = 0.01) were associated 
with increased risk of UC. 

The homozygous variant genotype of TLR9 1174 OA 
(rs352139) (OR unadj : 0.67, 95% CI: 0.47-0.95, p = 0.03) and the 
heterozygous genotype of LY96 -1625 OG (rsl 1465996) 
(OR unadj : 0.74, 95% CI: 0.57-0.96, p = 0.02) were associated 
with reduced risk of UC. In addition, both the homozygous and 
the heterozygous variant genotypes of TLR4 T>C (rsl 554973) 
(OR adj : 0.67, 95% CI: 0.48-0.94, p = 0.02), TMFA -308 OA 
(rsl800629) (OR unadj : 0.75, 95% CI: 0.58-0.98, p = 0.04), IL23R 
OA (rsl 1209026) (OR adj : 0.52, 95% CI: 0.29-0.94, p = 0.03) 
and PTPM22 1858 OA (rs2476601) (OR unadj : 0.71, 95% CI: 
0.52-0.96, p = 0.03) were associated with reduced risk of UC 
(Table SI and S2). No associations were found after Bonferroni 
correction for multiple testing. 

The variant allele of the polymorphisms have been shown to 
increase TMFRSF1A expression (TMFRSF1A -609 OT 
(rs4149570)), increase TLRs and TNF-ot mRNA levels {PPARG 
OG (rsl801282)), increase TLR2 levels (TLR2 OT 
(rsl816702)), increase MFKBIA expression (MFKBIA 2758 A>G 
(rs696)), increase MD-2 (LY96) and TNF-a levels {LT96 -1625 
OG (rsl 1465996)), increases TMFA expression (TMFA -308 O 
A (rsl 800629)), decrease IL-17 serum levels (IL23R OA 
(rsl 1209026)) and decrease TNF-ot serum levels (PTPM22 1858 
G>A (rs2476601)), respectively (Table 2). The biological function 
of the polymorphisms TLR4 T>C (rsl 2377632), TLR4 T>C 
(rsl554973), TLR9 1174 OA (rs352139) and IL10 OT 
(rs3024505) are unknown. 

Thus, polymorphisms associated with increased TMFRSF1A 
expression (TMFRSF1A -609 OT (rs4149570)), increased TLRs 
and TNF-ot mRNA levels [PPARG OG (rsl801282)) and 
increased TLR2 levels (TIR2 OT (rsl816702)) were associated 
with increased risk of UC. Furthermore, polymorphisms associat- 
ed with decreased IL-17 serum levels (IL23R G>A (rsl 1209026)) 
and decreased TNF-ot serum levels (PTPM22 1858 OA 
(rs2476601)) were associated with reduced risk of UC. In contrary, 
lower activity of the NFkB pathway through increased MFKBIA 
(iKBot) expression [MFKBIA 2758 A>G (rs696)), an inhibitor of the 



NFkB pathway, was associated with increased risk of UC. In 
addition, polymorphisms associated with increased MD-2 (LT96 — 
1625 OG (rsl 1465996)) and TNF-ot levels (TMFA -308 G>A 
(rsl 800629)) were associated with reduced risk of UC. 

Polymorphisms associated with risk of IBD 

The studied polymorphisms all showed the same direction of 
effect for both diseases except for the polymorphism in PPARG 
(rs 180 1282) (Table SI and S2). 

When including all patients (IBD), the homozygous variant 
genotype of TLR2 A>T (rs4696480) (OR unadj : 1.33, 95% CI: 
1.01-1.74, p = 0.04) and IL6R OT (rs4537545) (OR adj : 1.46, 
95% CI: 1.02-2.08, p = 0.04) and both the homozygous and the 
heterozygous variant genotypes of TLR2 OT (rsl816702) 
(OR adj : 1.42, 95% CI: 1.08-1.88, p = 0.01), TLR4 T>C 
(rsl2377632) (OR adj : 1.28, 95% CI: 1.01-1.64, p = 0.05), TLR9 
-1486 T>C (rsl87084) (OR adj : 1.29, 95% CI: 1.00-1.66, 
p = 0.05), TMFRSF1A -609 OT (rs4149570) (OR adj : 1.32, 
95% CI: 1.03-1.68, p = 0.03) and IL10 OT (rs3024505) 
(OR unadj : 1.25, 95% CI: 1.02-1.52, p = 0.03) were associated 
with increased risk of IBD. The homozygous variant genotype of 
TLR4 T>C (rsl554973) (OR unadj : 0.67, 95% CI: 0.46-0.98, 
p = 0.04) and both the homozygous and the heterozygous variant 
genotypes of IL2 '3 R G>A (rsl 1209026) (OR adj : 0.43, 95% CI: 
0.28-0.67, p = 2*10" 4 ) and PTPM22 1858 G>A (rs2476601) 
(OR adj : 0.62, 95% CI: 0.46-0.85, p = 3*10" 3 ) were associated with 
reduced risk of IBD (Table SI and S2). 

After Bonferroni correction for multiple testing both the 
homozygous and the heterozygous variant genotypes of IL23R 
OA (rsl 1209026) (OR adj : 0.43, 95% CI: 0.28-0.67, p = 0.007) 
and PTPM22 1858 OA (rs2476601) (OR unadj : 0.61, 95% CI: 
0.48-0.77, p = 0.001) were associated with reduced risk of IBD. 

The variant allele of the polymorphisms have been shown to 
increase IL-6r and IL-6 levels (IL6R OT (rs4537545)), increase 
TLR2 levels (TLR2 OT (rsl816702)), increase TMFRSF1A 
expression (TMFRSF1A -609 OT (rs4149570)), decrease IL-17 
serum levels (IL23R OA (rsl 1209026)) and decrease TNF-a 
serum levels (PTPM22 1858 G>A (rs2476601)), respectively 
(Table 2). The biological function of the polymorphisms TLR2 
A>T (rs4696480), TIR4 T>C (rsl2377632), TLR9 - I486 T>C 
(rsl87084), TLR4 T>C (rsl554973) and ILK) OT (rs3024505) 
are unknown. 

Thus, polymorphisms associated with higher IL-6r and IL-6 
levels (IL6R OT (rs4537545)), higher TLR2 levels (TLR2 OT 
(rsl816702)) and increased TMFRSF1A expression (TMFRSF1A - 
609 OT (rs4149570)) were associated with increased risk of IBD. 
In addition, polymorphisms associated with decreased IL-17 
serum levels (IL23R OA (rsl 1209026)) and decreased TNF-a 
serum levels (PTPM22 1858 OA (rs2476601)) were associated 
with reduced risk of IBD. 

Haplotype analysis 

Haplotype analyses of TLR2, TLR4 and TLR9 among patients 
with CD, UC and IBD versus healthy controls are shown in Table 
S3, S4, S5, respectively. Four haplotypes in TLR2, three in TLR4 
and two in TLR9 described 88%, 94% and 97% of the observed 
genotypes, respectively. 

The TLR2 haplotype combination 33 encompassing all the 
wildtype alleles (rs4696480AA, rsl 1938228CC, rsl816702CC and 
rs3804099TT) was associated with reduced risk of CD (OR: 0.19, 
95% CI: 0.08-0.45, p = 0.00007) and combined CD and UC (OR: 
0.32, 95% CI: 0.17-0.60, p = 0.0005). 

No associations were found for TLR4. 
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Table 2. The biologic effect of the studied single nucleotide polymorphism (SNP) and odds ratios (OR) for polymorphisms which 
have been shown to be associated with risk of Crohn's disease (CD), ulcerative colitis (UC) or inflammatory bowel disease (IBD) in 
previous studies and in this study. 



Previously found 
associations. Disease, 
genotype, OR (95% CI), 

Gene (SNP) rs-number Effect of the SNP p-value 



Associations found in this 
study. Disease, genotype, 
OR (95% CI), p-value 



TLR2 (activates inflammation through the canonical NFkB pathway) 


-16934 A>T rs4696480 Unknown [14] 


ND 


IBD: TT, 1.33 (1.01-1.74), 
p = 0.04 A 


OA rs11938228 Unknown [14] 


ND 


No association 


OT rs1816702 rsl 81 6702T increase 

receptor level c [31] 


ND 


CD: TT, 2.36 (1.08-5.16), 
p = 0.03 A ' B ;UC:CTorTT, 1.46 
(1.10-1.93), p = 0.009 AB ; IBD: 
CT or TT, 1.32 (1.05-1.65), 
p = 0.02 AB 


597 T>C rs3804099 597C decrease TNF-a, IL-1 fj 

& IL-6 level E [32] 


ND 


No association 


TLR4 (activates inflammation through the canonical or non-canonical NFkB pathway) 


G>A rs5030728 Unknown [14] 


ND 


No association 


T>C rsl 554973 Unknown [14] 


ND 


UC: TC or CC, 0.67 (0.48-0.94), 
p = 0.02 B ; IBD: CC, 0.67 
(0.46-0.98), p = 0.04 A 


T>C rsl 2377632 Unknown [14] 


ND 


UC: TC or CC, 1.42 (1.01-2.00), 
p = 0.04 B ; IBD: TC or CC, 1.28 
(1.01-1.64), p = 0.05 B 


TLR5 (activates inflammation through the canonical NFkB pathway) 


1 1 74 OT rs57441 68 1 1 74T (392 STOP ), decrease 

TNF-a, IL-1 B & IL-6 level G [32] 
and inhibit TLR5 function D ' E [33] 


CD: CT, 0.14 (0.03-0.57), 

p = 0.002 (Jewish) [28]; CD: No 

association (Non-Jewish) [28,29] 


No association 


TLR9 (activates inflammation through the canonical NFkB pathway) 


-1486 T>C rs187084 -1486C&1 174G decrease 

expression 0 [34] 


ND 


IBD: TC or CC, 1.29 
(1.00-1.66), p = 0.05 B 


1174 G>A rs352139 -1486C&1 174G decrease 

expression 0 [34] 


ND 


UC: AA, 0.67 (0.47-0.95), 
p = 0.03 A 


LY96 (MD-2 binds to and is involved in the TLR2 or the TLR4 complexes) 


-1625 OG rs11465996 -1625G increase MD-2 & 

TNF-a leve1 0E [35] 


ND 


UC: CG, 0.74 (0.57-0.96), 
p = 0.02 A 


CDU (binds LPS and transport it to TLR4) 


-159 G>A rs2569190 -159AA increase CD14 

level E [36,37] 


IBD: GA or AA, 2.95(1.77- 
4.90),p = 2*10~ 5 [22] (Korean) 


No association 


MAP3K14 (NIK is a central kinase in the non-canonical NFkB pathway) 


T>C rs7222094 rs7222094CC decrease 

NIK activity E [38] 


ND 


No association 


SUM04 (SUM04 conjugates to kBa and negatively regulates NFkB transcriptional activity) 


163 T>C rs237025 163C increase NFkBI 

expression 0 [39] 


ND 


No association 


NFKBIA (IkBa is an inhibitor of NFkBI) 


2758 G>A rs696 2758A increase expression 0 [40] 


UOExtensive colitis (Hungarian) 
[23]; CD: Inconclusive [24,25] 


UC: GA or AA, 1 .28 
(1.00-1.65), p = 0.05 A ' B ; CD: No 
association 


T>del rsl 7103265 rsl 71 03265del decrease 

expression 0 [41] 


ND 


No association 


NFKBI (NFkBI (p50/65) is a transcription factor. The NF-kB p50 subunit can act both pro-inflammatory as part of the p50/p65 com 
homodimer [42]) 


plex or anti-inflammatory as p50 


-94 ins/del rs28362491 -94del decrease p50 subunit 

expression 0,0 [43] 


Inconclusive [26] 


CD: Ins/- or -/- , 0.80 
(0.65-1.00), p = 0.05 A 


TNFA (TNF-a is a pro-inflammatory cytokine activated by NFkBI) 


-863 OA rs1800630 -863A increase 

expression 0 ' 5 [44] 


IBD: AA, 4.82 (2.60-8.96), 
p = 1«10~ 4 (Indian) [45] 


Failed to genotype 
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Table 2. Cont. 



Gene (SNP) rs-number Effect of the SNP 


Previously found 
associations. Disease, 
genotype, OR (95% CI), 
p-value 


Associations found in this 
study. Disease, genotype, 
OR (95% CI), p-value 


-857 OT rsl 799724 -857T increase 

TNF-a level aE - F [46] 


Inconclusive [27] 


Failed to genotype 


-308 G>A rs!800629 -308A increase expression c ' D [47] 


Inconclusive [27] 


UC: GA or AA, 0.75 
(0.58-0.98), p = 0.04 A 


-238 G>A rs361525 -238A decrease expression" [48] 


Inconclusive [27] 


No association 


TNFRSF1A (TNF receptor 1 (TNFR1) binds TNF-a and initiates a kinase cascade) 


-609 G>T rs4149570 -609T increase expression 1 " [49] 


ND 


CD: TT, 1.41 (1.02-1.94), 
p = 0.04 A ' B ; UC: TT, 1 .49 
(1.04-2.13), p = 0.03 AB ; IBD: 
GT or TT, 1.23 (1.02-1.50), 
p = 0.03 AB 


TNFA1P3 (TNF-a rapidly induced expression of TNFAIP3/A20 which inhibit NFkB activation and TNF-a mediated apoptosis) 


OG rs6927172 rs69271 72G increase expression 013 

[50] 


nd 


No association 


IL1B (pro-inflammatory cytokine activated by NFkBI) 


-3737 G>A rs4848306 -3737A decrease transcription F 

[51,52] 


ND 


No association 


-1464 G>C rs1143623 rs1 143623C decrease IL-1 b level E 

[52,53] 


ND 


No association 


-31 T>C rsl 143627 -3 1C decrease 

expression"' 0 [52-54] 


No association (Danish) [3] 


No association 


IL1RN (IL-1 RA binds to the IL-1 receptor and inhibit IL-1 p signaling) 


T>C rs4251961 rs4251 961C decrease IL-1 RA level E 

[55,56] 


ND 


No association 


IL4R (IL-4 receptor, IL-4 significantly inhibit IL-1 7 production) 


A>G (I50V) rsl 805010 rsl 80501 0G increase IL-1 7 

level" [57] 


ND 


No association 


IL6 (pro- and anti-inflammatory cytokine activated by NFkBI) 


-6331 T>C rs10499563 -6331 C decrease 

expression 015 [58] 


ND 


No association 


IL6R (binds IL-6 and initiates a kinase cascade) 


OT rs4537545 rs4537545TT increase IL-6r and IL-6 

level but not TNF-a, IL-1 RA and CRP 
level E [59] 


ND 


CD: TT, 1.73 (1.12-2.66), 
p = 0.01 B ; IBD: TT, 1.46 
(1.02-2.08), p = 0.04 B 


ILW (activated by NFkBI, capable of inhibiting synthesis of pro-inflammatory cytokines such as IFN-y and TNF-a) 


-592 OA rs1800872 -592A increase expression 0 [60] 


No association (Danish) [3] 


No association 


OT rs3024505 Unknown [3] 


CD: T-allele, 1.12 (1.07- 
1.17),p = 2*10~ 14 [3,16,61]; 
UC: T-allele, 1.25 (1.19-1.32), 
p = 6*10~ 17 [3,17] 


CD: No association; UC: CT or 
TT, 1.42 (1.10-1.82), 
p = 0.007 A ; IBD: CT or TT, 1.25 
(1.02-1.52), p = 0.03 A 


IL17A (activated by NFkBI, pro-inflammatory cytokine, potent mediator in delayed-type 
reactions, induces the production of IL-1J3, IL-6 and TNF-a) 


197G>A rs2275913 197A increase expression"' 13 [62] 


ND 


No association 


/L23fl (IL-23 receptor, IL-23 induce the production of IL-1 7 and IFN-y) 


G>A (R381 Q) rsl 1 209026 rsl 1 209026GG increase IL-1 7 serum 

level E [63] 


CD: G-allele, 2.66 (2.36-3.00), 
p = 1«1(T 64 [16]; UC: 
G-allele, 1.74 (1.57-1.92) , 
p = 5*10~ 2s [17] 


CD: GA or AA, 0.39 
(0.26-0.59), p = 1*10~ 5A ' B ; UC: 
GA or AA, 0.59 (0.39-0.90), 
p = 0.01 A ' B ; IBD: GA or AA, 0.47 
(0.34-0.66), p = 9*10~ 6A ' B 


IFNG (IFN-y is a pro- and anti-inflammatory cytokine activated by NFkBI) 


874 T>A rs2430561 874A decrease IFN-y level E [64] 


ND 


No association 


TGFB1 (TGF-fJI is a cytokine which can inhibit the secretion and activity of many other cytokines 


ncluding IFN-y and TNF-a) 




-509 OT rs1800469 -509T increase expression 0 ' 5 [65] 


ND 


No association 


PTPN22 (involved in several signaling pathways associated with the immune response) 
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Table 2. Cont. 



Previously found 



Gene (SNP) rs-number 


Effect of the SNP 


associations. Disease, 
genotype, OR (95% CI), 
p-value 


Associations found in this 
study. Disease, genotype, 
OR (95% CI), p-value 


1858 OA rs2476601 


1858A decrease TNF-a in serum 
level E ' F [66] 


CD: G-allele, 1.26 (1.17-1.37), 
p = 5*1CT 9 [16] 


CD: GA or AA, 0.54 
(0.41-0.72), p = 2*10~ 5A ' B ; 
UC:GA or AA, 0.71 (0.52-0.96), 
p = 0.03 A ; IBD: GA or AA, 0.61 
(0.48-0.77), p = 4*10~ 5A ' B 


PPARG (PPAFty is a transcription factor) 


OG (Pro12Ala) rs1801282 


rs1801282G decrease PPARy 
mRNA level, but upregulations 
MyD88 TLR4, TLR5, TLR9, P65 
and TNF-a mRNA levels E ' F [67] 


CD:GG, 0.33 (0.12-0.94), p = 0.03 
(Hungarian) [68]; CD: No 
association (Danish) [18]; 
UC: GG, 2.30 (1 .04-5.08), 
p = 0.04 (Danish) [18] 


CD: No association; UC: GG, 
2.12 (1.01-4.45), p = 0.05 A 


NLRP3 (NALP3 is involved in the inflammasome) 


OT rs4612666 


rs4612666T decrease 
expression 0 [69] 


ND 


No association 



A Crude (unadjusted). 

B Adjusted for age, gender and smoking status. 

c Function examined by flow cytometry. 

D Function examined by luciferase reporter assay. 

^function examined by enzyme-linked immunosorbent assay (ELISA). 

F Function examined by reverse transcriptase PCR (RT-PCR). 

G Function examined by electrophoretic mobility shift assay (EMSA). 

ND: not determined. 

doi:1 0.1 371 /journal.pone.009881 5.t002 



The TLR9 haplotype combination 11 (rsl87084CC and 
rs352139GG) was associated with increased risk of UC (OR: 
1.53, 95% CI: 1.03-2.28, p = 0.04). The TLR9 haplotype 
combination 1 2 was associated with increased risk of IBD (OR: 
1.35, 95% CI: 1.04-1.76, p = 0.03). Indirectly, these results 
support the analysis of the individual SNPs, where the variant 
allele of rs352139 (included in haplotype combination 22 which 
was used as reference) was fond to be associated with lowered risk 
of UC. 

Discussion 

In this study of severely ill patients, 16 functional polymor- 
phisms in 1 3 genes involved in regulation of inflammation were 
found to be associated with CD, UC or CD and UC combined 
(IBD) (Figure 1). As shown in Table 2, four known susceptibility 
loci for CD (PTPN22 (1858 OA) [16]), UC (ILK) (rs3024505 C> 
T) [3,17] and PPARG (rsl801282 OG) [18]) or CD and UC 
(IL23R (rsl 1209026 G>A) [16,17]) were replicated. Eleven new 
polymorphisms associated with risk of CD (IL6R (rs4537545 O 
T)), UC (TLR4 (rsl 554973 T>C and rsl 2377632 T>C), TLR9 
(1174 OA), LY96 (-1625 OG), NFKBIA (2758 A>G) and 
PTPN22 (1858 G>A)) or CD and UC (TLR2 (-16934 A>T and 
rsl816702 OT), TLR9 (-1486 T>C) and TNFRSF1A (-609 
G>T)) were identified. Other cohort studies of patients with CD 
and UC have found that other polymorphisms in the TLRs (toll 
like receptors) [19] and TNFRSFs (TNF receptors) [16,17,20,21] 
were associated with CD or UC. 

The biological interpretation indicates that a genetically 
determined higher activity of the inflammatory genes TLR2 
(rsl816702 OT), TJVFRSF1A (-609 OT) and IL6R (rs4537545 
OT) was associated with increased risk of CD and lower activity 
of the inflammatory genes NFKB1 (NFkB) (-94ins/del ATTG), 
IL23R (rsl 1209026 OA) and PTPN22 (1858 G>A) was 
associated with reduced risk of CD (Table 2). The picture was 



less clear for UC. A genetically determined higher activity of the 
inflammatory genes TNFRSF1A (-609 OT), PPARG (rsl801282 
C>G) and TLR2 (rsl816702 OT) was associated with increased 
risk of UC and lower activity of the inflammatory genes IL23R 
(rsl 1209026 OA) and PTPN22 (1858 OA) was associated with 
reduced risk of UC. In contrast, a genetically determined higher 
activity of the NFkB inhibitor IkBoc (NFKBIA 2758 A>G) and 
lower activity of TLR9 (haplotype 11 (-1486CC and 1174GG)) 
were associated with increased risk of UC. Furthermore, a 
genetically determined higher activity of LT96 (—1625 OG) 
and TNFA (—308 OA) was associated with reduced risk of UC. 
However, the risk of CD and UC seem to have shared 
mechanisms through the TLR2 (rsl 8 16702 OT), TNFRSF1A 
(-609 G>T), IL23R (rsl 1209026 OA) and PTPJV22 (1858 O 
A) polymorphisms confirming that the inflammatory pathways are 
involved in risk of both CD and UC [16,17]. 

This study was unable to confirm the associations between the 
variant aUele of ILK) (rs3024505 OT) and increased risk of CD 
[3,16] or the associations between the variant allele of CD14 (— 
159 G>A) and increased risk of CD, UC and IBD [22]. This study 
may be underpowered to detect an association in IL10 (rs3024505 
OT), as a GWA study with more than six thousands CD cases 
and fifteen thousands control found that the varaint allele was 
associated with increased risk of CD but with an odds ratio of 1 . 1 2 
[16]. However, IL10 has previously been found to be associated 
with risk of CD and UC in the Danish population [3]. The 
difference of association in GDI 4 (—159 G>A) may be due to 
genetic differences between the Korean and Danish population. 

The variant allele in NFKBIA (IkBoc) (2758 A>G) [23-25] and 
TNFA (—308 OA) were associated with higher and lower risk of 
UC, respectively, and the deletion polymorphism in NFKB1 
(NFkB) (— 94ins/del ATTG) was associated with lower risk of CD 
in our cohort study. However, there seems to be no consensus 
regarding these polymorphisms in other cohort studies [24-27]. 
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Figure 1. Sixteen functional single nucleotide polymorphisms (SNPs) in 13 genes involved in regulation of inflammation were 
found to be associated with susceptability of severe Crohn's disease (CD), ulcerative colitis (UC) or inflammatory bowel diseases 
(IBD). Eleven of the SNPs have not previously been reported as susceptability polymorphisms of CD, UC or IBD (TLR2 (rs4696480 and rs1816702), 
TLR4 (rs1554973 and rs12377632), TLR9 (rs187084 and rs352139), LY96 (rs1 1465996), NFKBIA (rs696), TNFRSF1A (rs4149570), IL6R (rs4537545) and 
PTPN22 (rs2476601)). 
doi:10.1371/journal.pone.0098815.g001 



TLR5stop (1174 OT) has been found to be associated with 
reduced risk of CD in a Jewish cohort [28] but not in a non-Jewish 
and German cohort [28,29]. No association with TLR5stop was 
found in this cohort study. 

No associations were found for IL1B (-31 T>C) [3], IL-10 (- 
592 C>A) [3] or PPARG (rsl801282 C>G) [18] in accordance 
with our previous Danish IBD cohort study using independent 
cases but the same control group [3]. The other polymorphisms 
studied in MAP3K14 (rs7222094 T>C), SUM04 (163 T>C), 
TNFAIP3 (rs6927172 C>G), IL1RN (rs4251961 T>C), IL4R 
(rsl805010 A>G), IL6 (-6331 T>C), IL17A (197G>A), IFJVG 
(874 T>A), TGFB1 (-509 OT) and NLRP3 (rs46 12666 OT) 
were not statistically associated with CD, UC or IBD in our 
cohort. None of the polymorphisms studied in these genes have 
been reported to be associated with CD or UC. We can not 
exclude that we did not find associations between these 
polymorphisms and risk of disease due to lack of power in this 
cohort. 

The results in this study should be interpreted with care. TLR2 
(rs4696480 A>T), TLR4 (rsl554973 T>C), TLR9 (1174 G>A) 
and TGFB1 (—509 OT) were not in Hardy- Weinberg equilib- 
rium among the healthy controls which is probable due to chance 
because of the number of polymorphisms analyzed. When 
corrected adequately for multiple testing they did not deviate 
from Hardy- Weinberg equilibrium. In the light of the obtained P- 
values and the number of statistical tests performed, we cannot 
exclude that some of our positive findings may be due to chance. If 
the results were corrected for multiple testing only the well known 



susceptibility polymorphisms in IL23R (rs 11209026 OA) and 
PTPN22 (1858 OA) were associated with reduced risk of both 
CD and IBD. We successfully tested 37 polymorphisms and, 
assuming a 5% acceptance level, two would be expected to be 
associated with susceptibility by pure chance. In this study 16 
polymorphisms were found to be associated with susceptibility and 
the found associations were biologically plausible. A major 
strength was that this clinically homogeneous and well-charac- 
terised cohort was rather large including 1035 patients with IBD 
and 795 healthy controls. All the patients were considered for anti- 
TNF treatment and were therefore considered to have a severe 
disease course. Genetic determinants may be expected to be strong 
among severely ill cases [30]. 

In conclusion, 16 functional SNPs in 13 genes involved in 
regulation of inflammation were found to be associated with 
susceptibility of severe CD, UC or IBD. Eleven of the SNPs have 
not previously been reported as susceptibility polymorphisms of 
CD, UC or IBD (Figure 1) although other polymorphisms in most 
of these genes have previously been associated with susceptibility 
of CD, UC or IBD. Our results suggest that genetically 
determined high inflammatory response was associated with 
increased risk of CD and the large number of polymorphisms in 
the TLRs associated with risk of CD or UC support that the host 
microbial composition, diet or environmental molecules in the gut 
are important factors driving the inflammatory response in 
genetically susceptible individuals. 
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